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ABSTRACT. 

We highlight a new "transverse" interference effect, arising from the coupling between 
electrons and holes, induced by a superconducting island in contact with a normal metal. 
As an example we compute the electrical conductance G of a T-shaped mesoscopic 
sample, formed by joining a horizontal bar of normal metal to a vertical leg of the same 
material. With a superconducting island located on the vertical leg and the current 
flowing horizontally, we examine the dependence of G on the distance L of the island 
from the horizontal bar. Of particular interest is the differential conductance G(E) at 
an applied voltage E = eV, which due to quantum interference between electrons and 
Andreev reflected holes, is predicted to oscillate with both L and E. For a system with 
a spherical Fermi surface, the period of oscillation with E is i^EpjkpL and with L is 
nEp/kpE, where Ep is the Fermi energy and kp the Fermi wavevector. 



PACS Numbers. 72.10.Bg, 73.40.Gk, 74.50. 
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1. Introduction. 

Mesoscopic superconductors, formed by adding one or more superconducting islands 
to a phase coherent normal host, represent a new class of quantum systems, whose reali- 
sation has been made possible by recent advances in nano-lithographic techniques[l,2,3]. 
The possibility of such structures has led many theorists to consider incorporating su- 
perconductivity into modern theories of phase coherent transport [4- 13] and several ex- 
perimental groups to carry out related investigations[14-18]. In this paper, we analyze 
a simple demonstration experiment, carried out originally to address the prediction of 
superconductivity induced conductance suppression[12,13]. The experiment [18] measures 
the electrical conductance G of a T-shaped sample of size less than the quasi-particle 
phase breaking length, formed by joining a horizontal bar of silver to a vertical leg of the 
same material. With a superconducting island located on the vertical leg, a distance L 
from the T-junction and the current flowing horizontally, the dependence of G on the dis- 
tance L and on the applied magnetic field was measured. As the applied field was lowered 
and the island allowed to become superconducting, macroscopic conductance changes of 
arbitrary sign were obtained, in agreement with the predictions of references [12, 13]. 

In what follows, we compute the conductance G(eV) = 31/ dV of such a T-shaped 
"trombone" as the superconducting island is slid along the vertical leg. To our knowledge, 
this quantity has not yet been measured for such structures, although the work of 
references [14-18] demonstrates that G(eV) is experimentally accessible. We note that 
at zero temperature, in the strict linear response limit of vanishing applied voltage, the 
conductance oscillates with L on the scale of the Fermi wavelength A^. This rapid 
variation arises from normal reflection at the interface with the island and is found in 
quantum modulated transistors[19,20,21], with no superconductivity. In contrast, we 
predict that the differential conductance at finite voltage V is sensitive to quantum 
interference between electrons and Andreev reflected holes and show that this oscillates 
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with L on the much longer scale ~ rcE F /(EkF) : where E = eV , E F is the Fermi energy 
and kp the Fermi wavevector. Such an effect is reminiscent of Tomasch oscillations in 
tunnel junctions[22,23], although in the present context, the current flow is transverse to 
the direction of the geometric feature causing interference. 

2. Calculation of the electrical conductance. 

To demonstrate this new interference effect, one notes that if the ends of the 
horizontal bar of the trombone are connected to external reservoirs by normal, current 
carrying leads, then the two probe, differential electrical conductance is given by 

G{E) = e[dl/d(m - ^ 2 )]( Ml -^=i?), 

where /zi, [ii are the chemical potentials of the reservoirs. At zero temperature, the 
analysis of references [7,8], yields, in units of 2e 2 /h, 



Q{E) = 1 - R PP (E) + RhpjE) + a(T hh ,(E) - T ph ,{E)) 

1 + a 

where 

Rh P (E) + T h i p {E) 



a 



R P 'h'(E) + Tph'(E) 

and the notation of reference[7] has been adopted. The coefficients R PP (E) , Rh p (E) 
(Rh'h'(E), R p 'h'(E)) are probabilities for normal and Andreev reflection respectively 
for quasi-particles (quasi-holes) from the left (right) reservoir, while T p > p (E) : T h / p (E) 
(Tf l h'(E),T p f l r(E)) are the corresponding transmission probabilities. These are obtained 
by tracing over combinations of sub-matrices of the scattering matrix of the structure[7,8]. 

Consider first the T-shaped structure shown in figure 1, which comprises three 
one-dimensional wires, connected at a node. In the absence of superconductivity, the 
Bogoliubov - de Gennes equation 

H (r) A(r) \ / V(r) \ _ />(r) \ , . 

A*(r) -HZ(r_) U(r) ~ ^ U(r) i ^' 



reduces to two separate equations for particles and holes. For clean wires, in a region 
occupied by the jth wire, a solution of the particle (hole) equation consists of incoming and 
outgoing plane waves of amplitudes v k Aj, v k l ^ 2 Bj (v q 1 ^ 2 Aj, v q 1 ^ 2 Bj) respectively, 
where Vk (v q ) is the group velocity of a particle (hole) of energy E = h 2 k 2 /2m-E F = 
Ep — h 2 q 2 /2m. These satisfy equations of the form 



and 



bA /a 1 

B 2 = S(E) A 2 

bJ \a 3 



bA [a 1 

B 2 =S*(-E) A 2 
B 3 \A 3 



(3) 



(4). 



To investigate the possibility of geometric resonances, we now introduce a supercon- 
ducting boundary, a distance L along the vertical leg of the structure. For convenience, 
we consider an ideal boundary, with negligible normal scattering, which induces coupling 
between particles and holes only. If r a (E) (f a (E)) are Andreev reflection coefficients for 
particles (holes) incident on the interface, then the coupling takes the form 



A 3 = r a (E)B 3 A 3 =f a (E)B 3 



(5). 



This pair of equations allows us to eliminate the amplitudes A 3: A 3 from equations (3) 
and (4), to yield an equation of the form (see appendix) 



f B A 

B 2 

Bi 
\B 2 J 



Spp(E) SP h (E)\ 
S h P(E) S hh {E) ) 



( A A 

A 2 
M 
\A 2 J 



(6), 



where the 2x2 sub-matrices S i: >(E) satisfy the following particle-hole symmetry relations 
S hh (-E) = [S pp (E)]* and S ph (E) = —[S hp (—E)]*. All reflection and transmission 
coefficients are obtained from the modulus squared of elements of these sub-matrices. In 
particular, T p , p (E) = |Sff| 2 , R hp (E) = \S^\ 2 , T hh ,{E) = \S^\ 2 and R^h^E) = |Sg| 2 . 



As an example, consider a spatially symmetric scatterer, for which 

S(E) = t tj . (7) 

In this case, the conduction formula (1) reduces to G(E) = (R p i l (E)+R h i p/ (E)+T hh /(E) + 
T p f p (E))/2. and in the limit E/E F << 1, one obtains (see appendix) 

T plp (E)~\t\ 2 \l + r*t/d\ 2 (8), 

T hh ,(E)~\t\ 2 \l + rt*/d\ 2 (9), 

R hp ~Rp, h , ~ |t| 4 /|rf| 2 (10), 

where 

-|r| 2 -exp[-2i(fc-g)L] (11). 

In the above expressions, only the denominator d depends on the length L. At zero 

energy, since k — q~ k F E/E F: this dependence vanishes. At finite energy, there is a slow 
oscillation with period given by 

k F LE/E F = 7i (12). 



3. Numerical results in two-dimensions. 

The above analysis demonstrates that transverse geometric oscillations are a generic 
feature of structures of the kind shown in figure 1. For a more realistic interface, which 
induces both normal and Andreev reflection, it is clear from the analysis of references [19- 
21] that there would be an additional, more rapid oscillation on the scale of k F L. To 
illustrate this feature and to show that slow oscillations survive in more realistic systems 
of finite width, we now present the results of a numerical simulation of the 2-dimensional 
T-shaped structure, shown in figure 2 and described by a tight binding Bogoliubov - de 
Gennes operator of the form 

H=(%° _£.) (13). 
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In this equation H a nearest neighbour Anderson model on a square lattice, with off- 
diagonal hopping elements of magnitude unity and A a diagonal matrix with on site, 
particle-hole couplings. The latter are chosen to vanish in the normal regions and to have 
a value Ao in the region occupied by the superconducting island (shown black in figure 2). 
The scattering region is chosen to be M = 30 sites wide and 6 sites long and is connected 
to external leads of width M. Within the leads, the diagonal elements are equal to e 
and within the insulating portion of the scattering region (shown shaded) the diagonal 
elements are set to an arbitrarily large value e. Within the remaining (conducting) portion 
of the scattering region, diagonal elements of Hq are chosen to be random numbers, 
uniformly distributed between e — W and e + W. In what follows, the choice e = 0, 
e = 50, 7 = 1 was made and free end boundary conditions employed. In the experiments 
of reference [18], the conductance per channel was estimated to be of order 0.1 e 2 //i, which 
for the system of figure 2, is obtained with the choice W ~ 0.1. The results of figure 3 
were obtained with this choice of W and A = 0.2. 

As well as possible oscillatory dependencies on E and L, it is of interest to examine 
the effect of a magnetic field. Therefore for each set of random diagonal elements, the 
conductance G(E) was first obtained with zero magnetic field, then recomputed with a 
finite flux </> through the conducting region. For each choice of E and L, the average 
conductance G(E, L) was computed from results obtained for 100 different realisations 
of the random diagonal elements, with the same sets of random numbers being used for 
each E and L. The squares in figure 3 show results for = 0, while the circles show 
results for <fi = 6^o, where 4>q = hc/e. For the latter, the field was introduced through a 
Peierls substitution in the off-diagonal elements of H . 

At zero energy, the conductance oscillates on the scale of the lattice spacing, which 
for the system of figure 2 is of order the Fermi wavelength Ap. At this energy, particle-hole 
symmetry ensures that there is no interference between electrons and Andreev reflected 
holes from the normal-superconducting interface. At finite energies, the wavelengths 
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of particle and holes differ by an amount proportional to E and therefore quantum 
interference on a length scale proportional to E~ x can occur. This new oscillatory effect, 
is distinct from more familar Tomasch[22] and Macmillan-Rowell[23] oscillations, since 
the latter occur in tunnel junctions, while in the structure considered here, there is no 
insulating barrier. In addition, the effect reported here occurs for E < A as well as 
E > A . 

To obtain a crude estimate of the periods of oscillation shown in figure 3, one notes 
that in the absence of disorder and superconductivity, a vertical, normal lattice of width 
two sites possesses a dispersion relation for particles of the form E = — 27[cos/c y a+cos/c"a] 
and for holes, E = +2'j[cosq y a + cos/c"a], where k™a = mi/3, n = 1, 2. Choosing E = 0, 
yields values for the allowed wavevectors kp on the Fermi surface given by k F a = tc/3 or 
37r/3. Hence for small E, E ~ 7sin/c-Ja (k y — q y )a = ^{k x — q x )aV3/2. and the period of 
oscillation obtained from equation (12) is given by 

(L/a)(£/ 7 )~27r/V3 (14). 

While this is a somewhat crude estimate, one notes that for E/^y = 0.2, this yields L ~ 18, 
which is of the same order as the period shown in the bottom graph of figure 3. 

Figure 3 shows that oscillations occur both in zero and finite magnetic fields and as 
a consequence, the magneto-conductance can change sign, if either L or E are varied. 
The oscillations predicted by figure 3 are a finite fraction of the total conductance and 
therefore should be observable experimentally. In practice it is not possibly to slide the 
superconducting island of a quantum trombone, without growing a separate sample with 
an independent set of random impurities. Consequently oscillations on the scale of 
are difficult to observe and would be measured as random sample to sample fluctuations. 
In contrast the longer scale variation in the differential conductance should be accessible 
experimentally. Indeed, by fixing the length L and varying E, the oscillations reported 
above should be observable in a single sample. For L = 11, this behaviour is illustrated 
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in figure 4, which shows a slow variation with E on a scale proportional to L~ x . For 
L = 11, equation (14) yields a the period of E/^y ~ 0.3, compared with a value E/^y ~ 0.2 
obtained from figure 4. 

4. Discussion. 

We have demonstrated that a superconducting island in contact with a normal 
mesoscopic conductor, produces geometric resonances, even if the superconductor is far 
removed from the current path. Such effects should be observable in the mesoscopic 
structures described in reference [18]. The above results, which are based on the 
assumption of a rigid order parameter, do not explain a weak magnetic field feature 
found in [18]. Figures 3 and 4 also show that the magneto-conductance can change sign 
with varying E and L, but the magnitude of the effect is much smaller than that observed 
in reference [18]. Since the theory of reference[12] predicts that the electrical conductance 
of a hybrid superconducting structure is sensitive to small changes in the magnitude of 
the superconducting order parameter, a possible explanation of the measured small field 
feature is that weak proximity tails present in the silver are modified by the applied 
field. To avoid such tails, it would be of interest to repeat these experiments, but with a 
magnetic normal host such as Nickel [24] 
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Appendix. 



Combining equations (3) and (5) yields B 3 = £33^3 + S31A1 + S3 2 A 2 = A 3 /r a (E). 
Similarly combining equations (4) and (5) yields S3 = S33A3+S31A1+S32A2 = A3/f a (E), 
where the i, j'th element of Sij(E) has been denoted S^ and the 2, jth element of S*j(—E) 
denoted Sij. These yield expressions for A3, A3 of the form 

A3 = axAx + a 2 A 2 + /Mi + (3 2 A 2 

and 

A 3 = /Mi + (3 2 A 2 + u x A x + a 2 A 2 , 

where 

«i = S 33 S 31 /d, a 2 = S33S3 2 /d, 

"1 = S 33 S 31 /d, a 2 = S 33 S 32 /d, 

Pi = S 31 /(dr a (E)), fo = S 32 /(dr a (E)), 

h = S 31 /(df a (E)), p 2 = S 32 /(df a (E)), with d = {r a {E)f a {E))- 1 - S 33 S 33 . 
Substituting these back into equations (3) and (4) yields equation (6), with 

gpp — ( + S\ 3 ai S\ 2 + Si 3 a 2 \ 
£21 + S 23 a\ S 22 + S 23 a 2 / 



gph 



S13P1 S\3f3 2 

S 2 30i S 23 (3 2 



ghh _ ( + Sl 3 a\ S± 2 + 5'l3Q;2 \ 

S 2 i + S23011 S 22 + S 23 a 2 J 



s h P = 



Si 3 pi Si 3 (3 2 
S23P1 S 23 /3 2 



For a N-S interface at a point xq along the x-axis, Andreev's approximation for 
reflection amplitudes is r a (E) = expz[— 9 — (f)+(k — q)xo], f a (E) = expz[— 6 + (f)+ (k — q)xo], 
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where <j) is the phase of the superconducting order parameter and 9 = tan _1 [(|A| 2 — 
E 2 ) 1 / 2 /E] +ivQ(—E). In view of particle-hole symmetry r a (E) = —f*(—E), as expected. 
For a symmetric scatterer, with S-matrices of the form 

r t t\ /fit 
S(E) = | t r t , S(E) = S*{-E) = It f t 

t t r J \t t f 



one obtains 



T p 'p(E) = |t| 2 |l + ai| 2 , 
T hh ,{E) = |t| 2 |l + a 2 | 2 , 



R hp (E) = \i\ 2 \^ 2 



and 



In these expressions, 



and 



R p , h ,{E) = \t\'\p 2 \ z . 

cti = ijda.1 = rt/d 
h = t/{dr a )fh = t/(dr a ) 

d=(r a (E)f a (E))- 1 -\r\ 2 . 



For small E, provided t, r are slowly varying compared with (E / Ep)kpL, one may 
evaluate r, t at E = and write tat* and f ~ r*. This yields equations (8) to (11) of 
the main text. 
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Figure Captions. 

Figure 1. A quantum trombone formed from three one-dimensional wires connected 
at a node. This figure shows the arrangement of incoming and outgoing amplitudes for 
particles. 

Figure 2. A quantum trombone formed from a two dimensional, tight-binding lattice 
of sites. The superconducting island is coloured black and located on the vertical arm of 
the trombone. The shaded area represents an insulating portion of the scattering region, 
while white area shows the weakly disordered conducting region. The current flows from 
left to right along the horizonatal arm of the T-junction. The external leads and scatterer 
are of width M = 30 sites. The horizontal arm is 6 sites long. The horizontal arm and 
vertical leg are each 2 sites wide. 

Figure 3. Results for the ensemble averaged differential conductance G(E, L) as a 
function of the position L of the superconducting island, for E = (top graph), E = O.I7 
(middle graph) and E = O.27 (bottom graph). In each case A = 0.2. The squares show 
results for zero applied magnetic field and the circles for an applied flux of Q(f>o passing 
through the normal part of the sample, (ie. the white plus shaded regions of figure 2), 
where 4>o = hc/e. 

Figure 4. As for figure 3, except that G{E, L) is plotted as a function of E. In this 
figure A = 0.2 and L = 11. 
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